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The behavior of cosmological evolution is studied in the case when a phantom fluid, that con- 
tributes to the universe accelerated expansion, is introduced in an F{R) model. At the early stages 
of the universe history, the dark fluid is seen to give rise to a deceleration of its expansion. For t 
close to present time it works as an additional contribution to the effective cosmological constant 
and, later, it produces the transition to a phantom era, which could actually be taking place right 
] now in some regions of the cosmos, and might have observable consequences. For t close to the 

. Rip time, the universe becomes completely dominated by the dark fluid, whose equation of state 

' is phantom-like at that time. Our model, which is able to reproduce the dark energy period quite 

' precisely, may still be modified in such a way that the epoch dominated by an effective cosmological 

constant — produced by the F{R) term and by the dark fluid contribution — becomes significantly 
shorter, what happens when a matter term is included. The dark fluid with phantom behavior 
gives rise to a super-accelerated phase, as compared with the case where just the viable F{R) term 
contributes, ft is also seen explicitly that an F{R) theory can be constructed from a phantom model 
in a scalar-tensor theory, in which the scalar field does not behave as phantom (in the latter case the 
action for F{R) would be complex). Promising F{R) models which are able to cross the phantom 
divide in a convenient way are constructed explicitly. 
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I. INTRODUCTION 



Our first aim in this paper will be to build a reliable cosmological model by using, as starting point, modified 
. gravity theories of the family of the so-called F{R) theories which comprise the class of viable models, e.g., those 
' having an F{R) function such that the theory can pass all known local gravity tests (see [l| -[III)- In the next section, 
, we will investigate cosmological evolution as coming from F{R) gravity. We will consider, for F(R), some candidates 
C ' to produce inflation and cosmic acceleration in a unified fashion. In particular, we investigate in detail the behavior 
<N : of F{R) gravity as the contribution of a perfect fluid. As a crucial novelty, an additional matter fluid will be included 
which may play, as we shall see, quite an important cosmological role. In fact, it may decisively contribute to the 
two accelerated epochs of the Universe, what is to say that we will study a model where dark energy consists of two 
0^ , separate contributions. The possibilities to obtain precision cosmology are enhanced in this way. 

■ Two main cases will be discussed: F{R) cosmology with a constant equation of state (EoS) fluid, and F{R) 
^ , cosmology in the presence of a phantom fluid. In the last case, a couple of speciflc example will be worked through 

in detail, namely one with a phantom fluid with constant EoS and, as a second example, a fluid with dynamical EoS 
. , of the type proposed in Refs. [lB| and In this case a dark fluid is present which has an inhomogeneous EoS that 

■ can depend on the proper evolution of the Universe, what opens a number of very interesting possibilities. 

As is well known, F(R) gravity can be written in terms of a scalar fleld — quintessence or phantom like — by redefining 
the function F{R) with the use of a scalar field, and then performing a conformal transformation. It has been shown 
that, in general, for any given F{R) the corresponding scalar-tensor theory can, in principle, be obtained, although 
the solution is going to be very different from one case to another. Also attention will be paid to the reconstruction 
of F{R) gravity from a given scalar-tensor theory. It is known that the phantom case in scalar-tensor theory 
does not exist, in general, when starting from F{R) gravity. In fact, the conformal transformation becomes complex 
when the phantom barrier is crossed, and therefore the resulting F{R) function becomes complex. We will see that, 
to avoid this hindrance, a dark fluid can be used in order to produce the phantom behavior in such a way that the 
F{R) function reconstructed from the scalar-tensor theory continues to be real. We will prove, in an explicit manner, 
that an F{R) theory can indeed be constructed from a phantom model in a scalar-tensor theory, but where the scalar 
field does not behave as a phantom field (in which case the action for F{R) would be complex). Moreover, we will 
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explicitly show that very interesting and quite simple F{R) models crossing the phantom divide can be constructed. 



II. FRAMEWORK 



Our first aim in this paper will be to construct reliable cosmological models by using, as starting point, modified 
gravity theories of the family of the so-called F{R) theories which comprise the class of viable models, e.g., those 
having an F{R) function such that the theory can pass all known local gravity tests. We now consider the action 
corresponding to one of these theories which, aside from the gravity part, also contains a matter contribution, namely 
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x^[R + FiR) + Lm] . (1) 

Here = 8ttG, and L„i stands for the Lagrangian corresponding to matter of some kind. Note that the first term 
in ll]) is just the usual Hilbert-Einstein action, and the F(R) term can be considered, as will be shown below, as the 
dynamical part of some kind of perfect fluid (see [ll] and [2^), which may constitute an equivalent to the so-called 
dark fluids. The field equations corresponding to action ([1]) are obtained by variation of this action with respect to 
the metric tensor g^^, what yields 

R^.. - \9tJ.uR + R^uF{R) ~ ]^9t,uF{R) + .g,,.nF'(i?) - V^V,F'(i?) = k^T^, . (2) 

Here the primes denote derivatives with respect to R. We assume a flat FRW metric, then the Friedmann equations 
are obtained as the 00 and the ij components. They take the form 

3i?2 = K^p,n - \P{R) + 3(i? + H)F'{R) - 18F"{R){H^H + HH) , 

-'dH^ -2H = K^p„, + \f{R) - {3H^ + H)F'{R) - DF'iR) , (3) 

where H{t) — a/a. Note that both equations are written in such a way that the i^(_R)-terms are put on the matter 
side; thus we may define an energy density and a pressure density for these _F(i?)-terms, as follows 

PFiR) = "^FiR) + 3{H + H)F'{R) - 18F"{R){H^H + HH) , 

PFiB.) = If{R) - i'3H^ + H)F'{R) - DF'iR) ■ (4) 

Then, the Friedmann equations ([3]) take a simple form, with two fluids contributing to the scale factor dynamics. 
From the energy and pressure densities defined in ([4]), one can obtain the EoS for the dark fluid, defined in terms of 
the F{R) components. This is written as 

- (3i72 + H)F'{R) - DF'iR) 
wf(b.} - _ i^(^) _^ 3(-^ ^ H)F'iR) + DF'iR) - VoVOF'(i?) 

PFiR) = -PFiR) + 2HF'iR) - VoV"F'(i?) . (5) 

The EoS, Eq. defines a fluid that de pen ds on the Hubble parameters and its derivatives, so it may be considered 
as a fluid with inhomogcneous EoS (see [l6j or [13] )■ In absence of any kind of matter, the dynamics of the Universe 
are carried out by the F(i?)-component, which may be chosen so that it reproduces (or at least contributes) to the 
early inflation and late-time acceleration epochs. In order to avoid serious problems with known physics, one has to 
choose the FiR) function in order that the theory contains flat solutions and passes also the local gravity tests (see 
[2]-0])- To reproduce the whole history of the Universe, the following conditions on the F(i?) function have been 
proposed (see [3]-[3|): 

i ) Inflation occurs under one of the following conditions: 

lim FiR) = -A, (6) 

R,^oo 

or 



R 



hm FiR) = ai?" . (7) 



In the first situation the FiR) function behaves as an effective cosmological constant at early times, while 
the second condition yields accelerated expansion where the scale factor behaves as a(i) ~ t^"/^. 
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FIG. 1: The function F{R) as given by (|10p for n — 2. We see that at the current epoch {R ~ Ro), F{R) behaves as a 
cosmological constant while for R <x (inflation) its values grow as a power law. 



ii ) In order to reproduce late-time acceleration, we can impose on function F{R) a condition similar to the one 
above. In this case the Ricci scalar has a finite value, which is assumed to be the current one, so that the 
condition is expressed as 

F{Ro) = -2Ro F'{Ro) - . (8) 



Hence, under these circumstances, the F{R) term is able to reproduce the two different accelerated epochs of the 
universe history. An interesting example that satisfies these conditions has been proposed in Ref. 
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This model is studied in detail in Ref. Q, where it is proven that the corresponding Universe solution goes, in its 
evolution, through two different De Sitter points, being one of them stable and the other unstable and which can 
be identified as corresponding to the current accelerated era and to the inflationary epoch, respectively. Thus, the 
model ^ is able to reproduce both accelerated epochs when the free parameters are conveniently chosen — as was 
done in Ref. — in such way that the the F(R) theory in ([9]) produces two de Sitter epochs and grateful exit from the 
inflationary stage is achieved. In the same way, from a similar example given in Ref. |4| , we can consider the function 

^ RnaR^_ 

This function, represented in Fig. 1, gives rise to inflation at the early stages of the Universe, assuming the condition 
([7]) holds, while for the current epoch it behaves as an effective cosmological constant. This is explicitly seen in Fig. 1, 
where function (|10p is represented for the speciflc power n — 2. 

Function (jlOp leads, at the current epoch, to a perfect-fluid behavior with an EoS given by Pf(r) ^ ~Pf(R)- In the 
next section such kind of F{R) functions will be considered, with the inclusion of a matter Lagrangian in the action 
([1]), and the corresponding cosmological evolution will be studied. We will see that both inflation and the current 
acceleration can indeed be produced by the F{R) fluid, provided other components are allowed to contribute too. 



III. COSMOLOGICAL EVOLUTION FROM VIABLE F{R) GRAVITY WITH A FLUID 



In this section we will discuss cosmological evolution as coming from i^(_R)-gravity. We consider the function F{R) 
as given by (jlOp . Functions of this kind have been shown to yield viable models which comply with all known local 
gravity tests (see, e.g., 3]), and they are good candidates to produce inflation and cosmic acceleration in a uniflcd 
fashion. We will study in detail the behavior of F(i?)-gravity, by considering it as the contribution of a perfect fluid 
in the way already explained in the preceding section. As a crucial novelty, an additional matter fluid will be here 
incorporated, which may play, as we shall see, an important cosmological role. It may decisively contribute to the two 
accelerated epochs of the Universe, what is to say that we study a model where dark energy consists of two separate 
contributions. Some examples of phantom evolution will be then discussed, in which the F{R) contribution acts as a 
cosmological constant and the additional fluid behaves as a phantom held, what opens again interesting new venues. 
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A. F(R) cosmology with a constant EoS fluid pm ~ Wmpm 

We consider a Universe governed by action ([T]), where the F{R) function is given by p^ . The matter term is 
represented by a perfect fluid, whose equation of state is p,„ = WmPm (where Wm =cons.) In this case, by considering 
this F{R) term as a perfect fluid — with energy and pressure densities given by ^ — the Friedmann equations reduce 
to Eqs. ([3]). For simphcity, we wiU study the case where n = 2; then, the function F{R) of our specific model reads 
(the study can be extended to other values of n without much problem) 

, , R'^(aR^~(3) , , 

First of all, let us explain qualitatively what the aim of this model is. For simplicity, we neglect the contribution of 
matter, so the first Friedmann equation yields 

2_ i?2(ai?2-/3) ^ ^^^2 , i^^2R{a-iR'' - 2aR^ - li) 



ZB^ = — -1 -U- + 3(i/^ + B) 



2(H-7i?2) ' ^ ' > (l + 7i?2)2 



- X^F'^R)i^B-B + BB) + + ^(f3^ - - , (12) 

(1 + ^R^Y 

where R — 6{2B^ + B). It can be rewritten as a dynamical system (see Ref. Q): 

B^C, C^Fi{B,C), (13) 

and it can be shown that its critical points are those which give a constant Hubble rate {H — 0), i.e., the points that 
yield a de Sitter solution of the Friedmann equations. We can now investigate the existence of these points for the 
model (jlip by introducing the critical points B = Bo in the Friedmann equation (|12p 

= " 2(1 + 7i?2) + (l+7i?2)2 ■ (14) 

To simplify, it can be rewritten in terms of the Ricci scalar Rq = 12Bq, what yields 

^Rl - + IrI + ii?o = . (15) 

This can be solved, so that the viable de Sitter points (positive roots) are found explicitly. A very simple study of 
Eq. (fT5|) . by using Descartes' rule of signs, leads to the conclusion that Eq. (fT5|) can have either two or no positive 
roots. As shown below, one of these roots is identified as an effective cosmological constant that produces the current 
accelerated expansion of the Universe, while the second root can produce the inflationary epoch under some initial 
conditions. Then, the model described by the function ((TT|) is able to unify the expansion history of the Universe. 
In order to get a grateful exit from the inflationary epoch, stability in the vicinity of the critical points needs to be 
studied: the corresponding de Sitter point during inflation must be unstable. This can be achieved, as is already 
known for the function Q, by choosing specific values of the free parameters. Even in the case of stable dS inflation, 
exit from it can be achieved by coupling it with matter, by the effect of a small non-local term (or by some other 
mechanism) . 

Let us now study the details of this same model at the current epoch, when it is assumed that F(R) produces the 
cosmic acceleration and where the matter component is taken into account. Function (jlip is depicted in Fig. 1; its 
minimum is attained at R = Rq, which is assumed to be the current value of the Ricci scalar. Further, F(R) as given 
by pT|) behaves as a cosmological constant at present time. Imposing the condition /37/a >> 1 on the otherwise free 
parameters, the values of Rq and F{Rq) are then given by Q 



i?o - — , F\Rq) = , F{Rq) = -2Rq ^ . (16) 

For simplicity, we shall study the cosmological evolution around the present value of the Ricci scalar, R = Rq, where 
([TT|) can be expressed as F{R) ~ ~2Rq + fo{R~ Rq)^ + O [{R — RqY) , the solution for the Friedmann equations ([3]) 
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can be written as H{t) = Ho{t) + 5H, where at zero order the solution is the same as in the case of a cosmological 
constant, namely 

Hit) = i/^coth ( (1 + ^^»)V^, 
V V 

As pointed out in Rcf. ^4], the perturbations SH around the current point R = Rq may be neglected. Therefore, we 
can study the evolution of the energy density ^ for the F{R) term as the EoS parameter defined by ^ around the 
minimum of the F{R) function, which is assumed to be the present value of the Ricci scalar. For such purposes, it 
is useful to rewrite the Hubble function (fT7)) as a function of the redshift z, instead of t. The relation between both 
variables can be expressed as 




l + z 



a 
ao 



A sinh 




(18) 



where ag is taken as the current value of the scale parameter, and — PomO'Q^^^^^''"\ being pom the current value 
of the energy density of the matter contribution. Hence, the Hubble parameter pT|) is expressed as a function of the 
redshift 



(19) 



Thus, the model characterized by the action ^ with the function (jlip and a matter fluid with constant EoS, depends 
on the free parameters (a,/3,7) contained in the expression of F{R) and also on the value of the EoS parameter 
(wm)- When imposing the minimum value for the function F{R) to take place at present time (z = 0), the free 
parameters can be adjusted with the observable data. Then we study the behavior of our model close to z = 0, 
where the contributions of non-linear terms produced by the function pip are assumed not to modify the solution 
(|19p . In spite of the Hubble parameter being unmodified, for z close to zero, the energy density Pf{r) will in no way 
remain constant for small values of the redshift. To study the behavior of the F{R) energy density given by @ it is 
convenient to express it as the cosmological parameter, ^pi^n) = j§if^, which can be written as 



FjR) 
6iJ2(z) 



1 + 



Hjz) 



F'{R) - 18F"{R) 



H{z) 



Hjz) 
H{z) 



(20) 



where the Ricci scalar is given by i? = 6[2iJ^(z) + H{z)]. This expression (|20|) can be studied as a function of 
the redshift. As we are considering the solution p^ . which has been calculated near z = 0, where the F{R) 
function has a minimum, the second term in the expression (|20p is negligible as compared with the other two terms 
{F" {Rq)Ro / F' (Rq) ~ (3^ /a » 1). This aproximation, as the solution p9)) . is valid for values of the Ricci scalar 
close to i?o, where the higher derivatives of F{R) are small compared with the function. The aproximation is no 
longer valid when the F(R) derivatives are comparable with F(R). Then, as discussed above, the free parameters of 
the model can be fitted by the current observational values of the cosmological parameters, and by fixing the minimum 
of F{R) to occur for z ~ 0. 

We shall use the value for the Hubble parameter Hq — 100 /i km s~^ Mpc~^ with h = 0.71 ± 0.03 and the matter 
density 17^ = 0.27 ± 0.04 given in Ref. [2l[. In Fig. ^ the evolution of the F{R) energy density (^0]) is represented 
for the model described by (|lip . where the matter content is taken to be pressureless (wm — 0). The cosmological 
evolution shown corresponds to redshifts from z = 1.8 to z = 0. For redshifts larger than z = 1.8 perturbations around 
the solution ()17p are non-negligible, and the expression for the Hubble parameter (I17p is no more valid. However, in 
spite of the fact that the evolution shown in Fig. [5] is not a complete picture of the F{R) model given by pTjl . it still 
provides an illustrative example to compare it to the standard AC DM model, which is also represented in Fig. [2l 
around the present time. As shown in Fig. 2, both models have two common points for z = and z — 1.74, while the 
behavior of each model between such points is completely different from one another. This result shows the possible 
differences between F{R) models of the type pl}|) and the AC DM model, although probably other viable models 
for F{R) gravity may give a different adjustment to the ACDM model (see H-L^])- Furthermore, the effective EoS 
parameter for the F{R) fiuid, Wf(_r): given by the expression (0), is plotted in Fig. 3, again as a function of redshift. 
It is shown there that Wpi^ji) is close to —1 for z = 0, where the F{R) fluid behaves like an effective cosmological 
constant, while it grows for redshifts up to z = 1.5, where it reaches a zero value. According to the analysis of 
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observational dataset from Supernovae (see Ref.[22l), the results obtained for the evolution of the EoS parameter, 
represented in Fig. 3, are allowed by the observations. 

As a consequence, the F{R) model given by pTjl . and where the F{R) fluid behaves as an effective cosmological 
constant, is able indeed to reproduce the same behavior at present time as the ACDM model. On the other hand, as 
was pointed out in the section above, the F{R) model given by (fTTjl also reproduces the accelerated expansion of the 
inflation epoch, so that the next natural step to undertake with this kind of models should be to study their complete 
cosmological histories, and the explicit details allowing for a grateful exit from inflation, what should demonstrate 
their true potential. Also note that, in order to obtain a realistic well behaved model, further analysis should be carried 
out, as the comparison with the luminosity distance from Supernovae, or the data from CMB surveys, although this 
is a major task, that will be left for future work. Furthermore, F{R) functions of this kind might even lead to a 
solution of the cosmological constant problem, by involving a relaxation mechanism of the cosmological constant, as 
was proposed in Ref. [23]. The effective cosmological constant obtained could eventually adjust to the observable 
value. This will also require a deeper investigation. 
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B. F{R) cosmology in presence of a phantom fluid 



According to several analysis of observational data (see [25[ and |26|), the effective EoS parameter of the physical 
theory that governs our Universe should quite probably be less than —1, what means that we should be ready to cope 
with phantom behavior. This possibility has been explored in F{R) theory p^ . where the possibility to construct an 
F{R) function that reproduces this kind of behavior has been demonstrated, and where the possible future singularities 
envisaged, as the Big Rip, so common in phantom models, do take place. An example is given by (see [Toj ) 



/(i?)^i?W., ^here , ^ ^o + W^(^±TO)±T 



(21) 



here Hq is a positive constant, and f{R) = R + F{R). In this case, the solution, neglecting other possible components, 
is the following 

Hit) ^ ^ , (22) 

where the Hubble parameter diverges at the Rip time t^. Another kind of F{R) gravity that reproduces a phantom 
behavior, but avoids the Big Rip singularity, is 



f{R) ^[R- 6Ho{Ho + 1) + 5/2)] ^ ' ^ ^^^^ " , (23) 

where Hq is a positive constant (for other F{R) models with transient phantom behavior, see [ISl)- The solution 
obtained is 

H{t) = H,t+^^. (24) 

In this last example, the universe has a phantom behavior at large time but there is no Big Rip singularity. However, 
the main problem of those models is that the corresponding F{R) function is not well constrained at small scales. 

We will now study the behavior of the cosmological evolution when a phantom fluid is introduced that contributes 
to the accelerated expansion. A phantom fluid can be described, in an effective way, by an EoS Pph — WphPph, where 
Wph < ^1- This EoS can be achieved, for example, with a negative kinetic term or in the context of scalar-tensor 
theory. The microphysical study of this kind of fluids is a problem of fundamental physics that has been studied in 
several works (see for example [33|] and [34]), as it is also the study of stability of the solutions. But it is not the aim 
for this paper to discuss, in this depth, the features of phantom fields, and it will suffice for our purpose to consider 
the effective EoS that describes a phantom fiuid. The F{R) function we will consider is the same as above, given by 
PH)) . which will contribute as an effective cosmological constant at present time. The motivation for this case study 
comes, as will be further explored, from the difficulty one encounters in constructing viable F(R) functions which 
produce phantom epochs, and this problem has a translation in the scalar-tensor picture in the Einstein frame too. 
As recent observational data suggest (see [1^ and [11]), the effective EoS parameter for dark energy is around —1, so 
that the phantom case is allowed (and even still favored) by recent, accurate, and extensive observations. 

Example 1. First, we consider a phantom fluid with constant EoS, e.g., pph — WphPph, where Wph < —1 is a constant. 
As in the case above, the F{R) function will be given by (fTT|) . and it will be assumed that its current value is attained 

at the minimum i?o = (/3/a7)^^*. We here neglect other contributions, as the dust term studied previously. Friedmann 
equations take the form 

ii^' = yPp/. + ^ , H = -^Pph{l + Wph) . (25) 

An expanding solution (note that a contracting solution can be obtained from the above equations too) for these 
equations is given by 




H{t) - 7^ rr. + \hr, (26) 

2\l + Wph\[ts' 



where tg is called the Rip time, that means the instant where a future Big Rip singularity will take place. Although 
for t much bigger than the present time, the solution (|26p is not valid anymore — because perturbations, due to the 
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derivatives of the function F{R), become large — and as the Ricci scalar R = 6(2H^ + H) grows with time, this model 
will behave as F{R) ~ for large times, which is known to produce accelerated expansion, and whose behavior is 
described by (see 

^(*) « J^t ' ^^^^ 
where /iq — — , -, , . And the effective EoS parameter for large time is 

-e//=-l-^. (28) 

That is to say, the universe will go through a super-accelerated expansion stage due to the phantom fluid and to 
the contribution coming from F(R), until it reaches the Big Rip singularity. In this case where no other contribution, 
as dust matter or radiation, is taken into account, late-time acceleration comes from the phantom behavior when the 
dark fluid component dominates, while the Universe behaves as dark energy when the F{R) term is the dominant one 
(as was in absence of this phantom fluid). In other words, the universe would not accelerate as a phantom one alone, 
but it will as a dark energy fluid with Weff > — 1. This is a fundamental point: the additional dark fluid is essential 
for a universe described by the F{R) function given in to display a phantom transition. 

We can study the evolution of this phantom fluid by using the continuity equation: 

Pph + iHpph{l + Wph) ^0 , (29) 

which can be solved, at the present time, when the Hubble parameter is expressed by Eq. Then, the following 

solution is obtained 

Pph = Poph _ ^^9/^ , (30) 

where poph is an integration constant. As we can see, the energy density for the phantom fluid grows with time until 
the Rip value is reached where the energy density becomes infinite. On the other hand, the evolution of the F{R) 
term may be studied qualitatively by observing the expression given for the energy density of this fluid in , so that 
this evolution is similar to the one in the above case. That is, as the value of the Ricci scalar increases with time, it is 
natural to suppose that, in the past, the energy density belonging to F{R) had smaller values than at present, so that 
the matter dominated epoch could occur when the F{R) fluid and the phantom fluid were much less important than 
they are now. Again, for this kind of model the F(R) contribution amounts currently to an effective cosmological 
constant which drives the universe's acceleration. 

Example 2. As a second example of a phantom fluid, we consider one with a dynamical EoS of the type proposed in 
Refs. pjj] and Here, a dark fluid is present which has an inhomogeneous EoS that may depend on the proper 
evolution of the Universe. This kind of EoS can be derived from the dynamics of an scalar field with some characteristic 
potential and a variable kinetic term (see [13 )i or either it may be seen as the effective EoS corresponding to the 
addition of various components that fill up our Universe. An EoS of this kind can be written as 

Pph^wpph+9{H,H,H...;t) , (31) 

where w is a constant and g an arbitrary function. The interesting point is that it is possible to specify a function 
g so that a complete solution of the Friedmann equations is obtained. In this case, our aim is to study a fluid that 
at present (or in the near future) can behave as phantom; to that purpose we choose w = —I, and the role of the g 
function will be to determine when exactly the phantom barrier is crossed. Thus, our model will be described by the 
phantom fluid in (PT]) and the function F{R) of (dJ), while the matter component can be neglected. As an example, 
the EoS for the dark fluid is given by 

Pph = -Pph - [^h'{t) +Pf{r) + PFiR)^ , (32) 

where the prime denotes derivative with respect to time, and Pf{r) and Pf{r) are the energy densities defined in 
As it is shown, through the EoS defined in ((32)) . the dark fluid will contribute to the acceleration of the Universe 
as a dark energy, and subsequently as a phantom fluid when it crosses the barrier Wph < — 1. We thus see that in 
this model accelerated expansion comes from two contributions, the F{R) term PT|) and the dark fluid one ([5^ . so 
that the accelerated expansion stage could cross the phantom barrier if the dark fluid dominates and contributes as 
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a phantom one. The difference with respect to the former model is that in the present case the dark fluid changes 
its behavior in the course of the expansion history (this will be seen again in an example below). As the EoS given 
in ([32|) can be rewritten in terms of the Ricci scalar R, it can be seen as additional terms to our F(R) function, in 
order to get the transition to the phantom epoch in the context of F(R) gravity. With the EoS ([5^ . the Friedmann 
equations can be solved, the following solution being found 



H{t) = h{t) 



(33) 



Different solutions can be constructed by specifying the function h(t). We are here interested in those solutions that 
give rise to a phantom epoch; for those cases the following splitting of the Hubble parameter is relevant 



Hi 



(34) 



where Hq, Hi and ts (Rip time) are positive constants. This function describes a Universe that starts in a singularity, 
a.t t = — which may be identified with the Big Bang one — then evolves to a matter dominated epoch, after which 
an accelerated epoch starts which is dominated by an effective cosmological constant and, finally, the Universe enters 
into a phantom epoch that will end in a Big Rip singularity. As in the cases above, the condition (|16p remains, that 
is F'{Ro) = 0, where Rq is the present value of the Ricci scalar. Fig. 4 is an illustration of how this model works: 
we see there that the Universe goes through a decelerated epoch until it enters a region where the evolution of its 
expansion has constant Hubble parameter and the F{R) term behaves as an effective cosmological constant. Finally, 
at tph the Universe enters into a super-accelerated phase which is dominated by the dark fluid of Eq. ([5^ . until it 

eventually reaches the Big P- — „j- j. _ j. rr\.„ c j-i,:„ ™„ j„i :„ 4.t,„4. j-u„ „ — 3sing phantom barrier takes 

place very softly, as it is se( of the Universe expansion. 

Alternatively, the evolution |^ ,s 



dH(t)/dt 
F'(R) 



tph ts 

t 

FIG. 4: The first derivative of tlie Hubble parameter is shown here to illustrate the different epochs of the Universe evolution 
(full curve). Also, the first derivative of the F{R) term with respect to R is represented as a function of time (dotted curve); 
the constant Hubble parameter region is given by a constant F{R), which plays the role of an effective cosmological constant. 
Starting at t — tph , the Universe enters into a phantom region dominated by the dark fluid, which ends at t = ta when the Big 
Rip singularity takes over. 



Wph 



rh'jt) +Pf(r) + Pf(r) 
^h^{t)-PFiR) 



(35) 
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and we may look at its asymptotic behavior, 



2{Ho + 4) 

For ^<t«h~^..,n--l + ^^y 

4 ■ 

For t ^ tQ — > ^ -I . 



2 

Forio, iph « t<ts — > Wph 1 - ^:i^Hl-^H^+ls,\ ' ("^^^ 
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We thus see that, at the starting stages of the Universe, the dark fluid contributes to the deceleration of its expansion. 
For t close to the present time, t^, it works as a contribution to an effective cosmological constant and, after t = tph, 
it gives rise to the transition to a phantom era of the cosmos, which could actually be taking place nowadays at some 
regions of it. Our hope is that it could even be observable with specific measurements. Finally, for t close to the 
Rip time, the Universe becomes completely dominated by the dark fluid, whose EoS is phantom at that time. This 
model, which is able to accurately reproduce the dark energy period, may still be modified in such a way that the 
epoch which is dominated by the effective cosmological constant, produced by the F{R) term and by the dark fluid 
contribution, becomes significantly shorter. This is supposed to happen when a matter term is included. 

To finish this section, the inclusion of a dark fluid that behaves as a phantom one gives rise to a super-accelerated 
phase, as compared with the case where just the viable F{R) term contributes. In the two examples here studied, we 
have proven that, while the F{R) term contributes as an effective cosmological constant, the dark fluid contribution 
produces the crossing of the phantom barrier, and it continues to dominate until the end of the Universe in a Big Rip 
singularity. In other words, both the contribution of F(R) and of the phantom fluid are needed. This is very clearly 
seen before, and the nice thing is that, because of this interplay, we have shown the appearance of nice properties of 
our model that no purely phantom model could have. Specifically, for our models above, F{R) gravity together with 
phantom matter, the effective w value becomes in fact bigger than -1, so that we are able to show that F{R) gravity 
can solve the phantom problem simply by making the phantom field to appear as a normal one. 



IV. SCALAR-TENSOR THEORIES AND F{R) GRAVITY WITH A FLUID 

We now turn to the study of the solutions given above in the alternative, and more commonly used, scalar-tensor 
picture. In Refs. [ll], [l^], it was pointed out that F{R) gravity can be written in terms of a scalar 

field — quintessence or phantom like — by redefining the function F{R) with the use of a convenient scalar field and 
then performing a conformal transformation. The scalar-tensor theory thus obtained provides a solution which is 
characterized by this conformal transformation, whose expression depends on the precise form of the F{R) function. 
It has been shown that, in general, for any given F(R), the corresponding scalar-tensor theory can in principle be 
obtained, although the solution is going to be very different from case to case. Also attention has been paid to the 
reconstruction of F{R) gravity from a given scalar-tensor theory. It is also known (see [l^) that the phantom case 
in scalar-tensor theory does not allow, in general, the corresponding picture in F{R) gravity. In fact, the conformal 
transformation becomes complex when the phantom barrier is crossed, and then the resulting F{R) function becomes 
complex too. To avoid this hindrance, a dark fluid can be used, as in the models of the preceding section, in order to 
produce the phantom behavior in such a way that the reconstructed F{R) function continues to be real. This point 
is important, and will be clearly shown below. Also to be remarked is the fact that scalar-tensor theories, commonly 
used in cosmology to reproduce dark energy (see for example fsS]), provide cosmological solutions whose stability 
should be studied in order to demonstrate the validity of the solution found. This has been investigated, e.g., in Ref. 
[s^ l and requires a very deep and careful analysis. For the purpose of the current paper, we will here concentrate in 
the proof the existence of the solution in the scalar-tensor counterpart, and the corresponding stability study will be 
left for future work. 

We start with the construction of the scalar-tensor theory from F{R) gravity. The action ([1]) can be written as 

S = J d^x^ [P{^)R + Q(0) + L™] , (37) 

which is known as the Jordan frame. Here F{R) has been written in terms of a scalar held. To recover the action ([1]) 
in terms of F{R), the scalar field equation resulting from the variation of the action ([57)1 with respect to is used, 
which can be expressed as follows 



P'(0)i? + Q'(0)=O, 



(38) 
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where the primes denote derivatives with respect to (j>. Then, by solving equation ([38)) we get the relation between the 
scalar field 4> and the Ricci scalar, (f> = (f>{R). In this way, the original F{R) function and the action ([1]) are recovered: 

R + F{R) = P{(t){R))R + Q(0(i?)) . (39) 

Finally, the scalar-tensor picture is obtained by performing a conformal transformation on the action (|37p . The 
relation between both frames is given by 

<?£m. = n^g,.. , where = P{cb) , (40) 

where the subscript e stands for Einstein frame. A quintessence-like action results in the Einstein frame 

1 



Ri 



where 
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d^PW) 



Q{4>) 

P2(0) 



and a{(j)) = P(0) , 



(41) 



are the kinetic term, the scalar potential and the coupling function, respectively. Hence, by following the steps 
enumerated above, we can reconstruct the scalar-tensor theory described by the action (I4ip for a given F{R) gravity. 
By redefining the scalar field (j) = R, and after combining Eqs. (|38p and ([39]) . the form of the two functions P(0) and 
Q{(f>) are found 



P(0) = 1 + F'(^) , Q(^)=P'(0)0-P(</)) 



(42) 



Hence, for a given solution in the Jordan frame ()37|) . the solution in the corresponding quintessence/phantom scalar 
field scenario — i.e., in the Einstein frame (j4ip — is obtained by the conformal transformation (|40p. and it is given by 



aE{tE)^[l + F'{(j>{t))f^^a{t) where dts ^ [I + F' {(j>{t))f^^ dt 



(43) 



We will be here interested in the phantom case. With that purpose, we analyze the model described in the above 
section by the F{R) function (fTTj) and the dark fluid with EoS ((32|) . and whose solution is (f34|) . For simphcity, we 
restrict the reconstruction to the phantom epoch, when the solution can be written as H(t) ^ and F{R) ^ ^R 
Using g2l), the function P(0(i)) takes the form P{t) - 2^i?2 ^ 2aHi{H,+i ) i 



frame. Then, using 



the solution in the Einstein frame is found 



(t.-t) 



ts-t ' -y- 

as a function of time t in the Jordan 



tE 



l2aHi{Hi + 1) 
7 



ln(is - t) — > aE{tE) 



l 2aHi{Hi + 1) 
7 



■ exp 



2aHi{Hi + 1) 



tE 



(44) 



Through the relation between the time coordinates in both frames, we see that while for the Jordan frame there 
is a Big Rip singularity, at t = t^, this corresponds in the Einstein frame to i —^ oo, so that the singularity is 
avoided, and there is no phantom epoch there. By analyzing the scale parameter in the Einstein frame, we realize 
that it describes a de Sitter Universe, while in the Jordan frame the Universe was described by a phantom expansion. 
As a consequence, we have shown that a phantom Universe in F{R) gravity may be thoroughly reconstructed as a 
quintessence-like model, where the phantom behavior is lost completely. This has been achieved in a fairly simple 
example and constitutes an interesting result. 

Let us now explore the opposite way. In this case, a phantom scalar-tensor theory is given and it is F{R) grav- 
ity which is reconstructed. As was pointed out in Ref. p^ . when a phantom scalar field is introduced, then the 
corresponding F{R) function — which is reconstructed, close to the Big Rip singularity, by means of a conformal 
transformation that deletes the kinetic term for the scalar field — is in general complex. As a consequence, there is no 
correspondence in modified gravity when a phantom scalar produces a Big Rip singularity. However, in our case we 
will analyze a scalar-tensor theory which includes a phantom fiuid that is responsible for the phantom epoch and for 
the Big Rip singularity. In this situation a real F(R) gravity will be generically reconstructed, as we are going to see. 

The action that describes the scalar-tensor theory is 



Se ^ d Xy^-gE 



Re - ^Lj{(f>)d^(j)d^'(l) -U{(l)) + a{<l))LphE 



(45) 
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where a{4>) is a coupling function and LphE the Lagrangian for the phantom fluid in the Einstein frame. In our case, 
we consider a phantom fluid with constant EoS, PphE = WphPphEi with Wph < —1- The Friedmann equations in this 
frame are written as 



Hi 



+ a(0)pp/i£;(l + Wph, 



(46) 



To solve the above equations, it turns out to be very useful to redefine the scalar field as (/) — t^- Then, for a given 
solution HE{tE), the kinetic term for the scalar field can be written as follows 



■^{He + 3(1 + Wph)Hl) - (1 + Wph)V{c^) 



1 - Wph 

To reconstruct F(R) gravity, we perform a conformal transformation that deletes the kinetic term, namely 



ff^ii^B = ^^fJ^Liy , where = exp 



(47) 



(48) 



Note that for a phantom scalar field, that is defined by a negative kinetic term, the above conformal transformation 
would be complex, as remarked in Ref. [18] and as will be shown below. Thus, the reconstructed action would be 
complex, and no F{R) gravity could be recovered. By means of the above conformal transformation, action (j45|l is 
given by 



S = d x^—g 



^kJ dtp-^ Lj(el>)^ 



R-e[^^^^^''^^]viq^)+Lph 



(49) 



where Lph is the lagrangian for the phantom fiuid in the Jordan frame, whose energy-momentum tensor is related 
with the one in the Einstein frame by = O^T^^^, and where we have chosen a coupling function q:((/)) = ^, for 
simplicity. By varying now the action ((49|) with respect to (j), the scalar field equation is obtained 



R 



(50) 



which can be solved as = 4>{R): so that by rewriting action ((49|l . the F{R) gravity picture result: 



F{R) 



-R 



(51) 



Hence, for some coupling quintessence theory described by action (P5|) . it is indeed possible to obtain a real F{R) 
theory, by studying the system in the Jordan frame through the conformal transformation (|48p. 

To demonstrate this reconstruction explicitly, an example will now be given. As we are interested in the case of a 
phantom epoch close to the Big Rip singularity, we will start from a solution in the Einstein frame He ~ ^ , and 
a scalar potential given by V{(t)) ^ {ts — 0)", with n > 0. Then, the kinetic term (|47|) is written as 



u;(0) 



-(2 + 3(l + ii;p,0) 



1 



1 - Wph [ts - 0)2 

The solution in the Jordan frame is calculated by performing the conformal transformation 

l/(lTfe/2) 



{ts - tE) 



ait) 



l±k 



(52) 



(53) 



where k — 



8(l+3(l+ttip/t)) 



and Wph < — 1. Note that, in this case, we can construct two diff'erent solutions, and 

It is easy to see that the 



correspondingly two different F{R) models, depending on the sign selected in Eq. 
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Big Rip singularity is thereby transformed, depending on the case, into an initial singularity (+), or into an infinity 
singularity (-). By using (|50p and (|5T|) . the following function F{R) is recovered 



To summarize, we have here shown, in an explicit manner, that an F(R) theory can be actually constructed from a 
phantom model in a scalar-tensor theory where the scalar field does not behave as a phantom one (in that case the 
action for F{R) would become complex). Moreover, very interesting F{R) models crossing the phantom divide can 
also be constructed explicitly [23|- The above reconstruction procedure, where we have taken the F{R) function of 
(fTTj) . can be generalized to other types of modified gravity models (for recent reviews see ^3^] and ,31,] )■ 



We have seen in this paper that the F{R) model given by pT|) . and where the F{R) fluid behaves as an effective 
cosmological constant, is able to reproduce the same behavior, at present time, as the ACDM model. On the other 
hand, this model gives rise to the accelerated expansion of the inflation epoch too, so that the natural next step 
to undertake with those models should be to study the complete cosmological history, in particular (what is very 
important) the explicit details for a grateful exit from inflation, what would demonstrate their actual potential. What 
is more, F{R) functions of this kind might even lead to a solution of the cosmological constant problem, by involving 
a relaxation mechanism of the cosmological constant, as was indicated in [23j (see also the very recent paper |32j]). 
The effective cosmological constant thus obtained could eventually adjust to its precise observable value. This issue 
is central and deserves further investigation. 

We have studied the behavior of the cosmological evolution when a phantom fluid is introduced that contributes 
to the accelerated expansion of the universe. The F{R) function we have considered is the one given in (|10p . which 
contributes as an effective cosmological constant at present time. The motivation for this case study comes, as it will 
be further explored, from the difhculty one encounters to construct viable F{R) functions which produce phantom 
epochs, and this has a rep resentation in the scalar-tensor picture in the Einstein frame. As recent observational data 
suggest (see [1^ and |26j|), the effective EoS parameter for dark energy is around —1, so that the phantom case is 
allowed — and actually favored by recent, extensive observations. 

We thus have seen that, at the early stages of the universe history, the dark fluid contributes to the deceleration 
of its expansion. For t close to present time, Iq, it works as a contribution to an effective cosmological constant and, 
after t ~ tph, it gives rise to the transition to a phantom era of the universe, which could actually be taking place 
right now in some regions of it. Our hope is that it could be actually observable. Finally, for t close to the Rip time, 
the Universe becomes completely dominated by the dark fluid, whose EoS is phantom-like at that time. This model, 
which is able to reproduce the dark energy period quite precisely, may still be modified in such a way that the epoch 
dominated by the effective cosmological constant, produced by the F(R) term and by the dark fluid contribution, 
becomes significantly shorter. This is the case when a matter term is included. 

The inclusion of a dark fluid with phantom behavior gives rise to a super-accelerated phase, as compared with the 
case where just the viable F{R) term contributes. In the two examples investigated in the paper we have proven 
that, while the F{R) term contributes as an effective cosmological constant, the dark fluid term produces the crossing 
of the phantom barrier, and it continues to dominate until the end of the universe in a Big Rip singularity. It is 
for this reason that the contribution of F{R) and of the phantom fluid are both fundamental. This has been clearly 
explained in the paper, and the nice thing is that, because of this interplay, we have shown the appearance of very 
nice properties of our model that no purely phantom model could have. This eliminates in fact some of the problems 
traditionally associated with phantom models and makes this study specially interesting. In particular, for some of 
our models of F{R) gravity together with phantom matter, the effective w value becomes in fact bigger than -1, so 
that we were able to show that F{R) gravity can solve the phantom problem simply by making the phantom field to 
appear as a normal field. 

To summarize, we have here shown, in an explicit manner, that an F(R) theory can indeed be constructed from a 
phantom model in a scalar-tensor theory in which the scalar field does not behave as a phantom one (in the latter 
case the action for F{R) would be complex). Moreover, very promising F{R) models which cross the phantom divide 
can be constructed explicitly (see also [13]). The above reconstruction procedure, where we have taken (|lip for the 
F{R) function, can be generalized to other classes of modified gravity models. 
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